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Abstract. We consider the question of existence of periodic solutions (called 
breather solutions or discrete solitons) for the Discrete Nonlinear Schrodinger 
Equation with saturable and power nonlinearity. Theoretical and numerical re- 
sults are proved concerning the existence and nonexistence of periodic solutions 
by a variational approach and a fixed point argument. In the variational ap- 
proach we are restricted to DNLS lattices with Dirichlet boundary conditions. 
It is proved that there exists parameters (frequency or nonlinearity parame- 
ters) for which the corresponding minimizers satisfy explicit upper and lower 
bounds on the power. The numerical studies performed indicate that these 
bounds behave as thresholds for the existence of periodic solutions. The fixed 
point method considers the case of infinite lattices. Through this method, the 
existence of a threshold is proved in the case of saturable nonlinearity and an 
explicit theoretical estimate which is independent on the dimension is given. 
The numerical studies, testing the efficiency of the bounds derived by both 
methods, demonstrate that these thresholds are quite sharp estimates of a 
threshold value on the power needed for the the existence of a breather solu- 
tion. This it justified by the consideration of limiting cases with respect to the 
size of the nonlinearity parameters and nonlinearity exponents. 

1. Introduction. This work concerns the Discrete Nonlinear Schrodinger Equa- 
tion (DNLS) 

i^ + e(A^) n -/? J F(|V'„| 2 )V'n = 0, /3 eR, (1) 
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on a finite lattice supplemented with Dirichlet boundary conditions, and on infi- 
nite lattices (n = (m, n^, ■ ■ ■ , n^) G Z N ). We concentrate on two examples of 
nonlinearities, 

F(\z\ 2 ) = and F(|z| 2 ) = \z\^ , (2) 

the saturable and power nonlinearity respectively. 

Note that we use the word power in two different senses in this paper, in power 
nonlinearity as above, and for a conserved quantity of the system fl} , defined as 

PM = £l^l 2 - (3) 

n 

We present some theoretical and numerical results related to the existence of 
time periodic solutions, having the form 

W^e^V* (4) 

Substitution of the expression into ([T]) with the nonlinearities © , shows that 
4> = {4'n}nei. N i satisfies the system of algebraic equations 

^ n = -e(A d </>) n + (3F{\4> n \ 2 )4> n , /Jet. (5) 

Solutions given by the expression (j4|), are called stationary wave solutions. Localized 
solutions fulfilling \<f> n \ — > as |n| — > oo, are known as discrete solitons or breathenQ- 
The problem of existence and properties of nonlinear localized modes in DNLS 
lattices, has attracted considerable research interest [5J[TT]. For recent studies on 
the saturable DNLS or its cubic-quintic approximation, we refer to [3) [19l HU 
[U (T5J [5D]. In these references, as well as in [U [T5] for the continuum models, 
remarkable properties and differences between models with power nonlinearities are 
reported. Although the case of the fundamental localized solutions assumes that <j) n 
is real [3U] , the results that we present here consider the existence and nonexistence 
of nontrivial breather solutions where <p n is in general complex. The existence 
of nontrivial breather solutions for DNLS |T|) will be established by variational 
methods. More precisely, we apply direct variational methods 2J to appropriate 
constrained minimization problems. This approach has been used to the focusing 
./V-dimensional DNLS equation with a power nonlinearity 

iij n + e(A d V)n + /^r^n = 0, /3 > 0, (6) 

in infinite lattices [21) . The results of [H] not only establish the existence of non- 
trivial breather solutions, but in addition the existence of a global minimum - an 
excitation threshold - in one of the fundamental conserved energy quantities, the 
power (or norm). This minimum requires the nonlinearity exponent to be greater 
than or equal to a certain critical value, depending on the lattice dimension. More 
precisely, it is proved in J5Tj Theorem 3.1, pg. 678], that if < a < jf, spatially lo- 
calized solutions (j4|) with Q < of arbitrary small power exist, while if a > there 
exists a ground state excitation threshold "Pthrcsh- The result of [21], resolved the 



1 We wish to reserve the term "breathers" throughout the text, for localized solutions of the 
form although this term is not strictly valid in the case of a finite lattice: this is due to the 
fact that the lower bounds for the power of solutions (|4j derived in this work in the case of finite 
lattice, are also valid for breather solutions. 
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conjecture for ground state breathers of [6]. The second conserved energy quantity 
associated with ©, is the Hamiltonian 

HM = <-^4>)2-^— V |0„| 2ff+2 , /3>o. 

(7+1 

nez" 

A ground state is a minimizer of the variational problem 

mf{HM :V[4>]=R 2 } , 

where (•, -)2 stands for the £ 2 -scalar product. The existence of the excitation thresh- 
old was proved with the help of a delicate discrete interpolation inequality similar 
to the Gagliardo-Nirenberg inequality of the continuous case. It is proved in [3T] 
that the inequality 

n6Z N \n£Z N ) 

holds for a > -S>. The excitation threshold is related to the best constant of ([7]). 
The ground state solution has frequency lu* and power Pthresh- 

After the preliminary results of Section 2, in Section 3 we consider the DNLS 
equation with saturable nonlinearity. Following the results mentioned above, for 
the DNLS equation with power nonlinearity (J6j) , by the application of the varia- 
tional approach to the saturable DNLS, we derive both the existence of nontrivial 
breather solutions, as well some bounds on the power of the minimizers. The vari- 
ational study considers the saturable DNLS, supplemented with Dirichlet boundary 
conditions. Although this is a simpler case in comparison with the infinite lattice 
(where one has to deal with the lack of compactness [21]), this case is of importance 
especially for numerical simulations. Since the infinite lattice cannot be modelled 
numerically, numerical investigations normally consider finite lattices with Dirichlet 
or periodic boundary conditions. We note that the choice of boundary conditions 
only matters, if a localized pulse is moving and collides with a boundary. 

In our study of the DNLS equation with saturable nonlinearity, we distinguish 
between defocusing (/? < 0) and focusing (f3 > 0) nonlinearity. For the defocusing 
case (Section 3.1), we consider two variants of minimization problems: seeking for 
nontrivial breather solutions <j) n {t) — e luJt (j) n , of prescribed frequency to > 0, in the 
first variant we consider a minimization problem for the energy functional 

£„[<!>] :=e(-A d (M)a+w 

\n\<K 

subject to a constraint on the logarithmic part of the saturable Hamiltonian of the 
DNLS CO), defined by 

W[0]=e(-Ad&0) 2 +/3 Y log(l + l^| 2 )- 
\n\<K 

This proves the existence of a nontrivial minimizer <f> of £ u (linear energy) and the 
existence of (3 < as a Lagrange multiplier, satisfying — (3 > u > such that 
4>n{t) = e luJt <fin is a solution of the saturable DNLS |T]) with this j3 as a nonlinearity 
parameter. We note that, in contrast to the DNLS with power nonlinearity ©, 
the frequency of the breather is limited by the condition A := —j3 > u, due to the 
resonance with linear modes. Also in contrast with the power nonlinearity case, 
the parameter f3 cannot be scaled out. Due to this fact, the result is of interest, 



4 



J. CUEVAS, J. C. EILBECK AND N. I. KARACHALIOS 



justifying the minimization of the linear energy, and the existence of a parameter (3 
for which this minimum is attained. 

The second variant for the defocusing case considers, for given (3 = —A < 0, the 
constrained minimization problem for the Hamiltonian 

M{H[4>] : Vl<j>] = R 2 } , 

that is, we study the existence of the nontrivial breather solution as a ground state. 
This approach proves the existence of a nontrivial minimizer </>*, at least in the 
parameter regime A > 2eN , and the existence of a frequency uj > (as a Lagrange 
multiplier), satisfying A > cu, such that ?/> n (£) = e lu3t (\>* n is a stationary wave solu- 
tion. Moreover, it is proved in this parameter regime that there exist frequencies, 
such that the corresponding nontrivial breather solution satisfy an upper bound for 
the total power, depending on the parameters A, e, N. The first numerical study 
performed on this parameter regime, for the behaviour of breather power0, justifies 
the existence of a range of frequencies, for which the upper bound of the power of 
the corresponding breather solution is satisfied. 

Section 3.2, is devoted to the focusing saturable nonlinearity (3 > 0. For this 
case, the existence of a nontrivial breather solution ip n (t) = e~ int (f> n , is proved 
similarly to Section 3.1, by considering the constrained minimization problem for the 
Hamiltonian (a ground state). Through the application of the variational method, 
a simple relation involving the frequency fi, /? and the power V[<f>] — R 2 is derived, 
which in terms of V, provides a local lower bound on the power of the minimizer. 
The result actually states that there exists f2 > satisfying this lower bound. We 
would like at this point to distinguish between the term excitation threshold in the 
sense of [21] and lower bounds, used throughout the text: The variational methods 
used here, establish the existence of minimizers of the energy functionals and the 
existence of parameters (/? or f2) as Lagrange multipliers which are associated with 
these minimizers. Furthermore they provide theoretically some local (in the sense 
that they depend on Q and j3) lower bounds on the power of the minimizers, with 
respect to the variation of these parameters and not the existence of excitation 
thresholds on the power in the sense of 21]. To this end, the numerical studies 
performed, investigate the efficiency of the lower bounds, as well their possible 
behaviour as the parameters fl and (3 vary. As a result, the second numerical study 
performed for the lower bound on the power of the focusing saturable case, verifies 
that this bound is actually the smallest value of the power below which we should 
not expect the existence of a nontrivial breather solution for arbitrary given to and 
(3 and dimension N. We call such smallest values thresholds on the existence of 
periodic solutions. 

In the case of the focusing saturable nonlinearity /3 > 0, the numerical studies 
verify that for some parameter values, the bound predicts the trend of the behaviour 
of the numerically computed power. Moreover the numerical study in 2D-lattices, 
shows that the breather solution of the focusing saturable DNLS demonstrates a 
similar behaviour to that of the focusing DNLS with power nonlinearity, with respect 
to the existence of excitation thresholds: power decreases as the frequency increases 
until it reaches a minimum value at a certain frequency, an excitation threshold in 
the sense of [21) . This behaviour should occur in higher dimensional lattices, and is 

2 We consider, in all the numerical studies throughout the paper, the power of single-site 
breathers (also known as Sievers— Takeno modes 18 ), which are the breathers with the smallest 
power 
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observed in [20] . It is worth taking into account that the saturable nonlinearity can 
be approximated by a power one with er = 1 for small values of \<j> n \- This relation 
holds when is close to (3 in the saturable case, and can justify the similarities 
between the saturable and power nonlinearities related to thresholds. 

In Section 4, we apply an alternative method to derive a threshold on the power 
of the breather solution of the saturable DNLS in the focusing case (3 > 0. We use a 
fixed point argument which was also used in [lOj . This approach is for the saturable 
DNLS, considered in infinite lattices (although similar estimates can be obtained 
in the case of Dirichlet boundary conditions). Replacing the saturable nonlinearity 
in the equation by its exact Taylor polynomial of order to, it is possible to derive 
the threshold for the power for both the DNLS with the saturable nonlinearity and 
the cubic-quintic approximation. The threshold appears to be the positive root of 
a polynomial equation. In the case of the ID-lattice, the numerical study verifies 
first that the power decreases as frequency increases, as it was predicted by the 
lower bound derived by the variational method of Section 3. Also, the numerical 
power approaches the predicted threshold derived by the fixed point argument as 
the frequency increases up to the limit (3. In comparison with the bound derived 
by the variational method, the latter also reaches the threshold derived by the 
fixed point argument, as the frequency increases. Especially for large values of the 
parameter f3, the theoretical estimates are proved to be quite sharp for large values 
of frequencies. For 2D-lattices we also observe numerically the appearance of the 
excitation threshold in the sense of [21]. In the focusing case, it follows from the 
numerical study that an increase of the dimension as well as of the nonlinearity 
parameter is required for this excitation threshold to appear. 

Section 5 is devoted to some theoretical and numerical results, related to the 
DNLS equation with power nonlinearity. We consider first the case of the defocusing 
(f3 > 0) DNLS, seeking for breather solutions ip n (t) = e~ iUt , > 0. By applying the 
same variational approach as for the saturable DNLS, we derive lower bounds on the 
power of the minimizers, depending on the dimension of the lattice. The numerical 
study demonstrates that these lower bounds can serve actually as thresholds on 
the existence of breather solutions. The approach of minimizing the linear energy 
functional 

£ n [<t>] :=e(-A d 0,0) 2 -ft £ |^„| 2 , Q > 0, 

\n\<K 

appears to be useful also in the power case as the numerical study indicates, since 
the derived lower bound on the power gives a result slightly closer to the real power. 

We would like to summarize by pointing out the main differences with the results 
of [21]. The results of [21] prove the existence of an excitation threshold on the 
power of periodic solutions of the focusing DNLS with power nonlinearity which 
appears for the case a > 2 /N, as well the existence of a frequency lo* > on 
which this threshold value on the power is achieved. The corresponding solution 
ipn{t) = e lu) <f> n is a ground state having power 7-thresh ~~ the excitation threshold 
value. 

The thresholds proved by the fixed-point argument are explicitly given threshold 
values which are independent of the dimension. There are explicit estimates satisfied 
for any periodic solution with frequency (3 > il > and for any N > 1 in the 
focusing saturable nonlinearity, and for any lu > in the case of the focusing power 
nonlinearity and for any a > 0, N > 1. They are thresholds in the sense that 
no periodic localized solution can have power less than the prescribed estimates. A 
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characteristic example for the justification of this claim as well as its usefulness, is 
provided by the numerical study regarding the focusing power nonlinearity. Even 
in the case a < 2/N, where the excitation threshold |21| does not exist, a periodic 
solution cannot have power less than the derived estimate. This "global character" 
of the estimates is revealed when one considers "limiting" cases of small (large) 
values of a < 2/N (large values of a > 2/N - the case of excitation threshold). The 
numerical studies verify that for small (large) values of frequencies the threshold 
is not only satisfied but is also a quite sharp estimate of the real power of the 
corresponding periodic solutions. Especially in the case a > 2/N, the numerical 
studies demonstrate the fact that the excitation threshold "Pthrosh, which is not 
known explicitly, satisfies the derived lower bound. Thus this bound should not be 
viewed as a prediction of the excitation threshold in the case of a > 2/N, nor as a 
theoretical prediction of the numerical power of periodic solutions, but as prediction 
of the smallest power a periodic solution can have, for any u,o~ and N > 1. The 
same global property is shared by the estimate for the saturable nonlinearity, as 
the numerical study in the case of large values of /3 in the defocusing case shows. 
We conclude with the remark that the numerical studies for the focusing saturable 
nonlinearity, as well as for the defocusing power nonlinearity, show that the local 
estimates derived by the variational methods also predict the smallest value of the 
power for arbitrary given values of parameters /3, a, uj and N. 



2. Preliminaries. For the convenience of the reader, we recall some basic results 
on sequence spaces and their finite dimensional subspaces as well on discrete oper- 
ators, that will be used in the sequel (see also pTU]. [21] ) . 

For some positive integer N, we consider the complex sequence spaces 

p>_\ § = {<M™ez«, n = {ni,n 2 , . . . ,n N ) £ Z N , </>„ € C, 

I M\ P = {J2n & » \<Pn\ P y <™ 

(The following elementary embedding relation [T7] holds between £ p spaces 

PcP, \\4>\\ P <M\ q 1<<?<P<oo, (9) 

in contrast with the L p (f2)-spaces, if f2 C M. N has finite measure. For p = 2, we 
get the usual Hilbert space of square-summable sequences, which becomes a real 
Hilbert space if endowed with the real scalar product 

(&^)a = Re <Mh, <t>, 4> e t 2 - (10) 

Note that any <f> n S i v , 1 < p < oo satisfies limijji—koo 4> n = 0, as assumed for spa- 
tially localized solutions (i.e. discrete solitons or breathers). The discrete Laplacian 
is defined for <j> n = <£(„ 1>rl2 ,..., njv ) as 

V 7 (ni,ri2 — l,...,njv) ~t~ ' ' " ~i~ Y(ni,Ti2r" 1) 
~^^' t P(n u n2.,...,n N ) + Y(m +l,n 2 . . ,n N ) 

(ii) 
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Now we consider the discrete operator V + : £ 2 — > £ 2 defined by 

(V + -0)™gZ« = {^(n 1 +l,n 2 ,...,n N )-i>(n 1 ,n 2 ,-,n N )} 
+ {V , (ni,n 2 + l,...,njv) — ^(ni ,712, ■■■,««) } 

- ?/>( rii ,ra 2 , ...,njv) }. (12) 

and V- : ^ 2 -> ^ 2 denned by 

(V~V>) neZ iV = {^(ni-l.na.-.njv) _ V , (m,r l2 ,...,njv)} 
+ {^(ni, n 2 -l, ...,njv) ~ V'(ni,n2,...,njv) } 

+ {V- 111 ,"2, } • (13) 

Setting 

ni,n 2 ,...,rajv )> (I 4 ) 
ni,ra 2 ,...,niv) ) 

(15) 

we observe that the operator — satisfies the relations 

N 

(-A d Vi,^ 2 ) 2 = ^(V+^,V+V- 2 ) 2 , for all e£ 2 , (16) 

(V+Vi>2) 2 = (^i,V-^) 2) for all Vi,^2 e^ 2 . (17) 

From (|T6|) . it is clear that — : I 2 — > £ 2 defines a self adjoint operator on £ 2 , and 
-A d > 0. 

To formulate the DNLS equation, subject to Dirichlet boundary conditions, we 
consider the finite dimensional subspaces of £ p for a positive integer K, defined by 

F(Z$) ={<f>et? : <t> n = for \n\ > K} . (18) 

We have £ P (Z^) = <C^ 2K + 1 ) N endowed with the norms ©-finite sums. In this case, 
for any 1 < p < q < oo, there exist constants C±, C% depending on K, such that 

Cl|M|p<||Vllg<C2|M|p. (i9) 

In the finite dimensional setting, the operator — A^ satisfies relations (|16p -(fTT |) . and 
its principal eigenvalue Ai > can be characterized as 

Al = inf M gJM . ,20) 

(I>e£ 2 {Z%) (M)2 ct>e£ 2 {Z^) 2-,|n|<jrl^»r 

Hence ([20]) and (12] ) -([13] ) imply the inequality 

N 

e\i £ |<^| 2 <e£||V+<^<4dV £ |^ n | 2 . (21) 

|n|<if f=l \n\<K 

Then it follows from that 

Ai < 4JV. (22) 



s 
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For example, in the case of an II? lattice n — l,...,K, the eigenvalues of the 
discrete Dirichlet problem — A<2</> — A</> with (j) real, are given by 



A„ = 4 sin^ — r , n = 1, . . . , K. 

For a N-dimensional problem, the eigenvalues are: 

A(ni,n 2 ,...,njv) 

2 ( Ti\H \ . 2 / n 2 7T \ . 2 / Tl N TT 

n ( v i(ZTT)J +sm 1 V 47ZTT)J + --- + Sm [w+i) ( 

for rij = 1, . . . , K, j = 1, . . . , N. In consequence, the principal eigenvalue of the 
discrete Dirichlet problem — A<2</> — A</> with <f> real, is given by 

Ai ee A(i,i ! ... > i) = 4ATsin 2 



3. Saturable nonlinearity: Constrained minimization problems-Dirichlet 
boundary conditions. 

3.1. A. Defocusing case (3 < 0: Periodic solutions i/j n (t) — e luJt (f) n ,ui > 0. In 

this section we consider the existence of breather solutions of the saturable DNLS 
equation, for the case (3 < 0. For convenience we set 

= -A, A > 0. 

Thus, we seek breather solutions for the DNLS equation 

i^ n + e(A d V0n+A ^" =0, A>0, (23) 

of the form 

ip n (t)=e iu)t <t> n ,u>0. (24) 
In this case, the system §5§ is rewritten as 

-e(A rf <ft) w +^ w =A - ■% A>0. (25) 

We note that in the case of the anticontinuum limit e = 0, it follows that the 
frequency of a non-trivial breather solution should satisfy 

A > lu. (26) 

The Hamiltonian TL and the power V, given by 

N 

n^] = e£||V+(fl|l-A E log(l + l^| 2 ), (27) 

v=l \n\<K 

m = E i<^i 2 > ( 2g ) 

\n\<K 

are quantities which are independent of time. We shall prove the existence of non- 
trivial breather solutions (|24p . by considering a constrained minimization problem. 
The system (|25|) will be considered as the Euler-Lagrange equation of the functional 

N 

^[0]=eEH V ^ll2+^ E l<^| 2 - A E l°g(l + l<^| 2 )> (29) 

v=\ \n\<K \n\<K 
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involving TL and V . To produce the Euler-Lagrange equation (|25|) from the func- 
tional TC and V, we shall use the following 

Lemma 3.1. Let <f> e £ 2 (Z^). Then the functional 

V{<f>) = log(l + |0«| 2 ), 

\n\<K 

IS a G 1 (£ 2 (Zf),R) functional and 

<V'(<^)=2Re E nn^^' ^ 2 (Z£). (30) 



Proof: We assume that (f>, i/j £ £ 2 (Z^). Then for any < s < 1, we get 
V((f> + sV) - V(0) _ 1 D „ y, f 1 d 

%\<K' 



-Re V / ^log(l + |«£ n + 0# n | 2 )d0 (31) 
2Re / / t , ijl — ; — 7T/ — prVrid-O. 



From the inequality 

E i p i^i^ E (I^i + W)i^i<(ihi^ + iiv»ii 2 )imi 2 , (32) 

we may let s — > 0, to get the existence of the Gateaux derivative (|3U)l (discrete 
dominated convergence) . 

To check that the functional V" : ^ 2 — > ^ 2 is continuous, we consider a sequence 
0m G ^ 2 (2<if) such that m — > in £ 2 (Z^). Then, by using an inequality very 
similar to ([93)1 (see Section [4]), we may show that \(V'((j) m ) — V'(4>), — > 0, as 
m — > oo. The result is also valid in the case of an infinite lattice (n £ Z w ). o 
The constrained minimization problem A.I. The first variational problem we 
shall discuss is a constrained minimization problem for the energy quantity 

AT 

SM := e Ell V ^H2+ w E l^l 2 ' ^W)' w>0 > ( 33 ) 

i/=l \n\<K 

for given w > 0. We have the following 

Theorem 3.2. Let u> > be given. Consider the variational problem on i 2 {JL^) 

inf | £ w [$ : E l0 S(l + \K\ 2 ) = R > \ ■ (34) 
[ \n\<K J 

There exists a minimizer <j> £ £ 2 (Z^) for the variational problem {3$ and a A(i?) > 
such that A(R) > uj, both satisfying the Euler-Lagrange eguation 

-e(A d 4>) n +u4> n = A - h \n\<K, 

l+\4>n\ 2 

<j> n = 0, \n\ > K. 
Moreover, it holds that 5Z| n |<jf log(l + |0n| 2 ) = R. 
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Proof: Clearly £ u defines a C 1 (£ 2 (Z^), R) functional. The minimization problem 
we intend to solve is for the functional £ u , restricted to the set 

Bi = I (ft e i 2 {K) : 53 Ml + I'M 2 ) =R>o\, (35) 

{ \n\<K J 

It is not hard to check that the sequence {</> m }meN £ B\ is bounded. Next, we 
consider a sequence {(ft m }meN € B±, such that <ft m — > (ft as m — » oo. We denote the 
nth coordinate of this sequence by (</> m )„. Using (fl"9)) . we observe that 



X) log(l + |(0 m )„| 2 )- ]T log(l + |0„| 2 ) 

M<K |n|<K 



<C||^ m -^|| £ i (zS) -0, (36) 



as TTi — > oo. Moreover, since <ft m € B\, we find from (|36|) that 

V log(l + \(ft n \ 2 ) = lim V log(l + \((ft m ) n \ 2 ) = i?. 



m — >oo 

|n|<i<" \n\<K 



Hence eft £ B\, which implies that B\ is closed. The functional £ UJ is bounded from 
below on £?i, since 



N 

|2 



£ u [cft] = e^HV^III.+c 1^1 

y=l |n|<K 

>W ^ |0„| 2 
|n|<K 

>w ^ log(l + |<£„ | 2 ) >w#. (37) 

|n|<K 

As we are restricted to the finite dimensional space £ 2 (Z^), it follows that any 
minimizing sequence associated with the variational problem (|34p is precompact: 
any minimizing sequence has a subsequence, converging to a minimizer. Thus £ u 
attains its infimum at a point (ft in B\ . We proceed in order to derive the variational 
equation for £ {JJ . Note that 

JV 

(^M,^> = 2 e 53(V+0,V+V) 2 + 2wRe 53 4>3 n - (38) 
y=l \n\<K 

By considering the ^(^(Zf ),K) (see Lemma ED]) 

£rW\ = 53 log(l + l0«| 2 )-i?, (39) 

|fi|<K 

we observe that for any (ft € Bi 

|n|<if 1 

The Regular Value Theorem ([2j Section 2.9], [9[ Appendix A,pg. 556 ]) implies 
that the set B\ = £^ 1 (0) is a C 1 -submanifold of ^ 2 (Z^). By applying the Lagrange 
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multiplier rule, we find the existence of a parameter A = A(i£) G M, such that 

JV 

I>=1 |n|<X 

-2ARc V — ^ — ^ n = 0, (41) 

\n\<K 1 + I"?"! 

for all -0 G ^ 2 (Z^). Setting ^ = in (JIT]), we find that 

N , ? ,2 

2^[0] = HV+0II1 + 2^ Yl \M 2 = 2A £ T^kj- (42) 

w=l |n|<K |n|<K 1 + 

Since <j> G -Bi cannot be identically zero and > 0, it follows from (|42|) that 

A > 0. 

Rewriting (142|) as 

e£||V+*||2 = A £ THX^-- E l^l 2 ' ( 43 ) 

!/=l \n\<K 1 + 10" I |n|<Jf 

we find that 

N 

e£l|V^<(A-^) J2 \4>n\ 2 - (44) 

The lhs of (|4"4"|) is positive, and we find that 

A(i?) > w. (45) 

Therefore condition ([2T?)) is justified. From (HT)) . there exists A > 0, such that the 
minimizer <fi S £>i solves the equation 

C-eAd^, V) 2 + ^(0, = A ( - , , forall VG^(Zf). (46) 



1 + |0| 



2 



The above formula, is clearly equivalent to the Euler-Lagrange equation ([231) . that 
is, any solution of (f2"S"| is a solution of PS)) , and vice versa, o 
The constrained minimization problem A. II. In the first minimization prob- 
lem, we derived, under sufficient conditions for given u> > 0, both the existence of 
(3 = —A < as a Lagrange multiplier, and the existence of a nontrivial <j> G ^ 2 (Z^), 
such that ([24)) is a breather solution of ([23)) . Here, by studying a different mini- 
mization problem, we shall derive some sufficient conditions depending on the given 
(3 = —A < 0, the lattice spacing e and the dimension of the lattice N. These pro- 
vide both the existence of a parameter u> > and a nontrivial tp* , such that ([55)) 
is a solution of (Q]) involving this uj as the frequency of the breather solution. This 
alternative variational approach for the existence of breather solutions ([24)1 for the 
DNLS equation ([2"3")h is to minimize the Hamiltonian 7i, constrained to the set 



B = { £ l\TL N K ) : £ l^l 2 = R2 > )• ( 47 ) 

\n\<K 
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Theorem 3.3. Let A, e, N > be chosen such that 

A > 4eN. (48) 
Consider the variational problem on £ 2 {lJ^) 

inf {H[4>] : V[<j>\ = R 2 > 0} . (49) 

Assuming that 

, A-4eiV , . 

* < ieiV ' < 50 > 

t/iere exists a minimizer <f>* G £ 2 (Z^) for the variational problem and w = 
> such that A > u>(R), both satisfying the Euler- Lagrange eguation 

C = 0, |n|>tf. 
Moreover, it holds that 7-* [</>*] = R 2 . 

Proof: We note first that H : £ 2 -> R, is bounded fr om below, since 

H[1>] >-A ^ log(l + |0„| 2 ) >-A J |0„| 2 = ~AR 2 . (51) 

|n|<if |«|<K 
Again, the finite dimensionality of the problem implies that any minimizing sequence 
associated with the variational problem (|4T[) is precompact, and that any minimizing 
sequence has a subsequence, which converges to a minimizer. Therefore, we conclude 
that Tt : B — > R attains its infimum at a point cf>* E B. The next step is to derive 
the variational equation ([5]). To this end, by using Lemma 13. 11 we observe that 

N 6 — 

(H'ttM) = 2e£(V^,V+V0 2 -2ARe ]T |2 V„- (52) 

v=\ \n\<K l<Pn| 

We consider next the functional Mr : ^ 2 (Z^) — > R 

Mr[4>] ■= ]T \^\ 2 -R 2 - (53) 

\n\<K 

Clearly, JVfl is ^(^(Zf ),R) and 

V) =2Re £ <^ n . (54) 

|n|<K 

Moreover, for any <f> E B, we have 

(A^M» = 2 J] I'M* = 2i?2 ^ °- (55) 

\n\<K 

Therefore, we may apply the Regular Value Theorem, to show that the set B = 
A^O) is a C^-submanifold of £ 2 (Z%). It follows from (55]), JM]) and the rule of 
Lagrange multipliers, that there exists a Lagrange multiplier £1 = fl(R) £ R, such 
that 

jv ,» 
(W'[^]-OA^[^]^> = 2e^(V+r,V+V0 a -2ARe £ 1 + |li2 ^ 

w=l |n|<K |(?)nl 

-20Re ^ #$^ = 0, for all tp e £ 2 '(Z$) . (56) 
\n\<K 
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2U(4>*) = 2n 2^ 10*1 = 2fti? 2 , (57) 



Thus, for t/j — (/>*, wc find that 

\cp | = zsi/ 32 

|n|<K 

where the functional W : £ 2 (Z^) — > M is defined by 

AT 



W(0) :=e£||V+<£||!-A E 



1+ 2 

For an appropriate choice of R, e, A, we can show that W[</>*] < 0: From (f2"Tj) . we 
have 

AT 

e£||V^|| 2 <4e7V|H| 2 . (58) 
1/=1 

Then, noting that |</>„| 2 < X)|n|<K 1 0^. | 2 , and using ([58]). we observe that 

Um < 4eN J2 - A E - (™ 2 - A^) ■ (59) 

\ti\<K \n\<K llVn.112 V / 

Therefore W[<£*] < 0, if 

4eiV(l + i? 2 ) < A. (60) 

which in terms of R, gives (|50l) . We shall consider equation ([57)1 . for the choice of 
parameters (1501) : since <fi* £ B cannot be identically zero, and U\(f>*\ < 0, it follows 
that Q(R) < 0. Thus, we may set 

Q(R) = -u(R), lu{R) > 0. (61) 

The inequality (|44|) is still applicable, to verify that uj{R) satisfies condition (|26|) . 
We get from ([TO|) and that there exists > 0, such that the minimizer 

cf>* e B solves the equation (|46f . o 

Numerical Study. The result of Theorem l3.2[ establishes the existence of a nontrivial 
0, for a given oj > 0, and the existence of A > 0, satisfying A > ui such that 
^n(£) = e luJt (/>n, solves the (l23l) equation, with such a A as a nonlinearity parameter. 
In Theorem 13.21 there are no further relations assumed between to, e, N. 

On the other hand, the result of Theorem 13.31 has the following implementation 
concerning the existence of breather solutions: for given A, e, N satisfying condition 
(BSD, 

A > 4eiV, 

there exists some ui > such that A > ui and <f>* G ^ 2 (Z^) not identically zero, 
solving the Euler-Lagrange equation ([23]) . Hence, if (|48p holds, there exists u> > 0, 
such that VVi(i) = e * w *^n) i s a solution for the DNLS equation (|23"|) , with power 

The statement of the Theorem 13.31 could be useful in the sense that the parameter 
regime (|48|) establishes the existence of a range for frequencies u> > 0, such that the 
corresponding breather solutions of the DNLS equation (f2"3"| have power satisfying 
the upper bound (f5TJ|) . 

A numerical study has been performed to study the behaviour of the power of 
the breather solution in the parameter regime (|48[> . to test the result of Theorem 
and the upper bound (|50p . There is an extra condition for the existence of 



14 J. CUEVAS, J. C. EILBECK AND N. I. KARACHALIOS 

Saturable nonlinearity. p=-2. N=1 . e=0.2 Saturable nonlinearity. p=-2. N=2. e=0.1 




A 



Figure 1. Power vs frequency of the defocusing DNLS with sat- 
urable nonlinearity for the cases (a) A = 2, e = 0.2, N = 1 (b) 
A = 2,e = 0.1, iV = 2. The predicted value by the upper bound 
(|50|) . in both cases is V = 1.5 (see text). In case (b), we observe 
the existence of a minimum of the total power (excitation thresh- 
old) 



breathers arising from the condition of non-resonance with linear modes, which is 
that 

A>4Ne-cu. (62) 

Clearly, (|62|) is satisfied when (|48|) is assumed. The numerical power verifies that 
there exists a range of frequencies such that the corresponding breather solutions 
have power satisfying the upper bound (|50[) : first we have depicted the power versus 
the frequency for breathers with A = 2 and e = 0.2 in a 1-dimensional lattice (see 
Fig. [J). From the theoretical prediction, an lu > should exist satisfying A > u>, 
with a power which should be always smaller than 1.5. From the figure, it can be 
deduced that the prediction is satisfied for all u> > 0.353. We have also considered 
the case A = 2, e = 0.1, for N — 2. Similarly, for the 2D-lattice, the breathers 
solution of frequency oj > 0.830, satisfy the theoretical upper bound V = 1.5. The 
numerical study in the 2D-case reveals the existence of an excitation threshold for 
the defocusing saturable DNLS as in the case of the power nonlinearity [21] : power 
decreases as frequency increases, attaining a minimum value for a certain value of 
frequency, as shown in the inset of Figure[T](b). As the frequency increases further, it 
seems that the power reaches a "threshold value" . In the case of higher dimensional 
lattices, the upper bound (|50p . could be even more useful, as an estimate from above 
of the excitation threshold as well as of the "threshold" value of the increased power 
reached, as the frequency increases further up to the resonant limit. 

3.2. B. Focusing case > 0: Periodic solutions ip n (t) = e" lf2t r>„,ft > 0. This 
section considers the existence of breather solutions of the saturable DNLS equation 
in the focusing case (3 > 0. We look for breather solutions of the form 



4> n (t) = e- mt (j) n , n > 



(63) 
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By considering again the case of the anticontinuum limit e = 0, it follows that the 
frequency of a nontrivial breather solution (163|) satisfies 

[3 > SI. (64) 

We choose to consider the minimization problem for the Hamiltonian 

N 

Hit] = eJ2Wt^\\i + E log(l + l<M 2 ), /3>0, (65) 

v=\ \n\<K 

constrained on the set B given by (|47p . since the approach of the variational problem 
A.I does not seem to be applicable in this case. 

Theorem 3.4. Let [3 > be given, and consider the following variational problem 
on£ 2 (Z%) 

M{H[<f>] : P[<j>] = R 2 > 0} . (66) 

There exists a minimizer cf> g £ 2 (Z^) for the variational problem f^ff| ] and SI = 
u)(R) > such that (3 > Sl(R), both satisfying the Euler- Lagrange equation 

-e(A^)„-0^ n = -0 h \n\<K, 

1 + \<Pn\ 

4> n = 0, \n\ > K. 
Moreover, it holds that V[(f>] = R 2 and 

n > iw ^ 

Proof: To see that Tt, is bounded from below, this time we use the inequality 
(|2"0f : we have 

N 

H[4>] > e£||V+0||l 

> eAi E l^| 2 >eA!i? 2 . (68) 

\n\<K 

The existence of the minimizer <f> € ^ 2 (^k )' an< ^ °f the Lagrange multiplier SI G K, 
can be derived by using the same variational arguments as in case A: since 

N 7 

H'[$]-n*r R [$M) = 2e^(V+^,V+^) 2 + 2/3Re ]T : Ti , 2 ^» 

f=l |ra|<K 1 + l^" I 

-20Re ^ 0„V^=O, for all r/> G ^ 2 (Z^). (69) 

|n|<if 

Setting ip = <fi, we get the equation 

N , 7 , 2 

2W[^]:=26X;i|V+0||3+2)9 £ TTiT^ = 20 £ M> ^ 

v=l \n\<K L + \ ( P n i \n\<K 

and since <f> ^ and > 0, we find that St(R) > 0. To justify condition we 
note by using ([58]) and (j70|) . that for arbitrary e > 

(4e7V + /3) £ |^| 2 >f7 £ |^„| 2 , 

\n\<K \n\<K 
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implying that 



8 0, 
€+ 4N > 4N' fOTany e> °' 



thus 

8>n>0. (71) 
Since V[4>] = R 2 , we get from (|7DJ}, that 



A' 



l&l 



|n|<K M<-R" 



< 20 ^ 



+ I0n| 2 



^ 1 + ||^ 112 



Let us assume that i? > and 



\n\<K \n\<K X ^ I IV" 1 12 

g 

1 + i? 2 

n<^. (73) 

Then, ([72]) and ([ST]) imply that i? = 0, and in this case ([73]) turns out that /3 < O. 
Thus, we have a contradiction both with the assumption R > as well as with ([51]) . 
Therefore for i? > 0, we should have /? > f2 and ([57]) . o 

On the other hand, the inequality (I67p . in terms of the power, it can be rewritten 

R 2 > £ - 1. (74) 



Thus the result of Theorem 13.41 shows that, for given 3 > 0, there exists some 
fl > 0, satisfying 3 > O, and a nontrivial € £ 2 (Z^), such that ip n (t) = e~ lQt (j> n , 
is a breather solution of fTJ, with power satisfying the lower bound (174]) . Let us 
observe that the rhs of (|74|) predicts that the power should be a decreasing function 
of the frequency f2, as the frequency increases to the resonant limit 3. In the next 
section, we shall derive a threshold value for the power of breather solutions for 
the focusing case 3 > 0, by using a fixed point argument. The lower bound of this 
section as well as the fixed point threshold, will be tested numerically. 

4. Thresholds for periodic solutions of the saturable DNLS by a fixed 
point argument-Infinite lattices: focusing case 8 > 0. We repeat here the 
fixed point argument of [10| to derive a threshold on the power, for the non-existence 
of non-trivial breather solutions for |T]). The approach covers the case of an infinite 
lattice (n £ TL N \ We consider the case where the parameters 8 > £1 > are given, 
and we investigate conditions on the non-existence of non-trivial solutions of the 
form 

^„W = e- i!i 'i, 8>n>0. (75) 
Note that 4> satisfies (0, rewritten as 

-e(A rf 0) ra -»0» =-P - -% 8>n>0. (76) 

1 + \<Pn\ 

For the convenience of the reader we state [33] Theorem 18. E, pg. 68] (Theorem of 
Lax and Milgram), which as for the case of the 2cr-power nonlinearity [10], will be 
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used to establish existence of solutions for an auxiliary linear system of algebraic 
equations related to (|7S|) . 

Theorem 4.1. Let X be a Hilbert space and A : X — > X be a linear continuous 
operator. Suppose that there exists c* > such that 

Re(Au,u)x > c*\\u\\x, for all u e X. (77) 

Then for given f £ X, the operator equation Au = /, u 6 X, /ias a unique solution 

Recall that for f(x) — j^x, x £ M, the following identity holds 



(_i \m+l J(m+1) 

/(x) = l-x 2 +a; 4 + --- + (-l) n a ; 2 " + ^ i = , for all x e R, (78) 

1 + ar 



which coincides with the Taylor polynomial of order m. Applying (|78p for x = 
\C\iC G C we rewrite the saturable nonlincarity as 

c 



no , , 

f_1\TO+l|/* 

I 2 J- l^l 4 J L Ln m l/*l 2m J- 1 I !i 

i + ICI 2 



i + ICP 

= (i-ici 2 + ici 4 + --- + (-inci 2m + v — iT!)i 

for all C £ C. Using (|79[) . equation ([75]) can be rewritten as 

- e(A d 0)„ + (f3 - CL)(t>n = F*{4> n ) + T.fa), (80) 

where 

F*(<t> n ) : = 0(\4> n \ 2 <Pn-\<Pn\Un + ---+(-l) m+1 \<t>n\ 2m ct>n), (81) 

T * (0n) = " TT^p • (82) 

Setting 

5:=f3-ft>0, (83) 
we observe that the (linear and continuous) operator Ts : i 2 — ► £ 2 , defined as 

(T s <t>) neZ N = -e(A d cf>) neId N + 6<j) n , (84) 

satisfies condition ([77)1 if (|83)) holds, since 

iv 

(T^» 2 =e^||V+0|| 2 + ^||0|| 2 > ( 5||0|| 2 for all <fi £ £ 2 . (85) 
i/=i 

Next, setting F 2m (C) = |C| 2m C, we may define from F 2m , a map F 2m : £ 2 — > I 2 . We 
have that 

||F 2m (z)||? 2 < J2 \z n \ 4m+2 <\\z\\i m+2 . (86) 

Then, writing F* = PJ^JLi^j, we may also define from a nonlinear map 
F* : I 2 ^ £ 2 , since 

m m 

||F.(*)||*» < W^MIU* < 0£ ( 87 ) 
i=i i=i 

Similarly, from the remainder term T*, we may define a map T* : £ 2 — > I 2 : we have 

</3 2 £ i^r +4 </3 2 ii0ii, 4 2 m+4 . 
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Hence the assumptions of Thcorcm l4.H are satisfied and the auxiliary linear problem 

[Ts4>)n& N = £( z n), (lC(z)) neZ N := (F*(2f)))„ eZ JV + (T* (z)) n &Z N ) ( 89 ) 

has a unique solution. 

We proceed with the definition of the map C on which the fixed point argument 
will be applied. For any given z £ £ 2 , we define the map £ : £ 2 — > £ 2 , by C(z) := <p, 
where cj) is the unique solution of the operator equation ([89| . Thus, C : £ 2 — > £ 2 is 
well defined. We will verify next that C : Br — > i3_R and is a contraction, satisfying 
the hypotheses of the Banach fixed point theorem. 

Let (, £ £ £\r such that <f> = £(£), tp = £(£). The difference \ := cj) — ip satisfies 
the equation 

(TsX)nez N = (^(C))«ez« - (£(0)nez N , (90) 

where 

(/C(C))„ e z« - (£(0)n 6 Z* = {(F,(C))n(EZ~ - (F.(0)n 6 Z"} 

+ {(T,(C))n 6 ^-(T,(0) n gz*}. 
The map F* : £ 2 — » ^ 2 is locally Lipschitz: We recall that for any F £ C(C, C) which 
takes the form F(z) = g(|C| 2 )C; with g real and sufficiently smooth, the following 
relation holds 



no - no 



{(C - 0(9(r) + rg'{r)) + (C- 0*V(0} dff, 



(91) 



for any C, £ G C,where $ = 6>C + (1 - 6% £ (0, 1) and r = |$| 2 (see [7, pg. 202]). 
Applying dlT]) for the case of F 2m {0 = |C| 2m C, one finds that 



F 2m (0 - F2m(0 = Km + 1)(C - 0l*| 2m + m(C - O$ 2 |*| 2m_2 ]d0. (92) 
Jo 

Assuming that C, £ € ^flj an d noting that ||$||2 < R,we get from the inequality 
J] |F 2m (Cn)-F 2m (e„)| 2 < /? 2 (2m+l) 2 J2 { f \^n\ 2m \Cn ~ Udo] 

< /3 2 (2m+i) 2 ^ \ f'wnincn-udo] 

< /3 2 (2m + l) 2 ^ ( f R 2m \C„ - £n\dd\ 



= /? 2 (2m + l) 2 R im 53 ic„-e„i 

nez™ 

Application of ([91]) to the remainder term, where 

( — l) m + 2 7' m + 1 



(93) 



1 + r 



implies that 

r,(c)-r„(0 = 
+ f\(~0® 2 



(C-0 



(-l) m+2 (m + 2)|$| 2m + 2 (_l)™+2|$|2m+4 



1+|$|2 (1 + |$| 2 ) 2 

(-l) m+2 (m + l)|$| 2m (_l)m+2|$|2m+2- 



l + l$l 2 



(1 + |$| 2 ) 2 



d0. 



d0 
(94) 
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Then working similarly as for the derivation of (f93|) . we get that 

|T4Cn)-T,(^)| 2 

2 

|2m+4 



<P 2 Y.{f t( 2m + 3)l$«| 2m+2 + 2|$„| 2m+4 J |C„ ~ UdO ! 

</3 2 E ( / [(2m + 3)||$||, 2 2 m+2 + 2||$||, 2 2 m+4 ] ICn-^l^l 

< /3 2 [[(2m + 3)i? 2m+2 + 2R 2m+i ] 2 £ |C„ - £„| 2 . (95) 



From inequalities (f93l) and (f95|) . we set 

m 

L\{R) = 5>j + l)i?^ 

i=i 

L 2 (#) = (2m + 3)i? 2(m+1) +2i? 2 < m+2 \ (96) 
i(i?) = Li{R) + L 2 {R). 

Then combining (|93[) and ([93]) , we observe that the map : i? 2 — > £ 2 is locally 
Lipschitz, satisfying 

\\K{Q-m)\\ P <pL{R)\\t-i\y. (97) 
Taking now the scalar product of (|9T)|) with \ in € 2 and using (|97p . we have that 

AT 

e£H V ^+ 5 IWl2 < ||/C(C)-/C(0ll2||xll2 



< (3L(R)\\C - Clbllxll 



< ^llxlli + S^WHC-^lli (98) 



From (f9"5)l . we obtain the inequality 



HxH 2 = ||£(z) - < ^L 2 (i?)||C - ^1 li- (99) 

Since £(0) — 0, we observe that the map £ : Br — > £>i? and is a contraction, having 
a unique fixed point - the trivial one - if 

L(R) < i (100) 

We consider the polynomial function 

Tl{R):=L{R)- 8 -. (101) 

A threshold value for the existence of nontrivial breather solutions can be derived 
from condition (I100|) : the polynomial equation H(R) = 0, has exactly two real roots 
R* < < such that i?» = -R*. Thus 

U(R) < for every R e (0,R*), 

that is, condition poop is satisfied if R £ (0,i?*). We summarize our results in the 
following 



20 



J. CUEVAS, J. C. EILBECK AND N. I. KARACHALIOS 



Theorem 4.2. We assume that the parameters e > and (3, Q > are given such 
that 

13 > n. 

Let > denote the unique positive root of the polynomial equation Tl{R) = 0, 
where n(i?) and L(R) are given by U01\) and i96\) respectively. Then a breather 
solution f75| ) of (QP, must have power bigger than 

V ,n ~ Rl(0,ty (102) 

Cubic-Quintic approximation. The saturable nonlinearity can be approximated 
by a cubic-quintic approximation (m = 2 in Taylor's formula (l79|) ). 

We consider first the approximation without taking into account the remainder 
term (the term L2(R) does not appear in the polynomial equation). Stationary 
wave solutions (175[) satisfy the infinite system of algebraic equations 

-e(A d( /») n -^, i = -/3(l-| ( /» n | 2 + | ( /» n | 4 )0„, neZ N , (3 > Q > 0. (103) 

In this case, the threshold value is Vp,n = i? 2 , where i?* is the root of the quadratic 

s_ 

a 



equation II(-R) = 3i? 2 + 5i? 4 - 4 = 0, S = [3 - SI, 



Setting, for example (3 = 2, fl — 0.5, we obtain the threshold value 

P 2 ,o.5 « 0.189898. 

Exact saturable nonlinearity. For the exact saturable nonlinearity we should 
take into account the remainder term: we look for the root R*, of the equation 
n(i?) = 3i? 2 + 5i? 4 + 7R 6 + 2R 8 - | = 0. For /? = 1, O = 0.5, we obtain 

^2,0.5 « 0.180917. 

Note that the threshold value appears to be the same for parameters (3 > > 0, 
giving the same ratio . 

Finite dimensional lattice. We may also derive a threshold value, taking into 
account the finite dimensionality of the lattice, when the problem is supplemented 
with Dirichlet boundary conditions. We may replace the constant S = (3 — > 
by the constant 

Si = eAi + /?-n, 

where (3 > Vt > 0. Therefore, in this case one has to work with the polynomial 
equation 

H(R) := L(R) - 5 -±. (105) 

P 

For the case of the cubic quintic approximation, the threshold value (|104p , (without 
considering the remainder term), becomes 

3 [29/3 + 20(eAi -O)] 1 / 2 
V W,n-~ + ^ , (106) 

involving the principal eigenvalue Ai > of the operator — and lattice spacing 
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Saturable nonlinearity. p=1 . N=1 Saturable nonlinearity. (5=1 . N=2 




Figure 2. Dependence of the power with respect to the frequency 
for ID and 2D lattices of the focusing DNLS with saturable non- 
linearity for (3—1 and e = 1. Power decreases approaching the 
theoretical lower bounds when Q, tends to (3. 

Numerical study. We performed numerical studies to test the lower bound (|74|) 
derived in Theorem OTHl and the threshold (|102|) derived in Theorem l4.2l In figures 
[Hand [31 we show the dependence of the power Vn = I 'An 1 2 with respect to the 
frequency ft. In the figures, the solid line represents the numerically calculated 
power, while the dashed line corresponds to the analytical threshold defined as the 
root of the equation (|105j) for the finite with Dirichlet b.c. lattice. The dot-dashed 
line corresponds to the lower bound defined in ([74]) . In most of the cases considered 
in our study, the difference between infinite and finite lattices is difficult to see, 
so the former has not been represented in the figures. As we cannot model an 
infinite lattice numerically, the numerical results for the power correspond to finite 
lattices with Dirichlet b.c. Note that, due to the resonance with linear modes, 
the frequencies of the breathers are limited by the condition f2 < j3. Besides, the 
continuation of the solutions is quite difficult close to this limit. 

Figure [2] refers to the parameters [3=1 and e = 1. The figure verifies that 
we should not expect existence of breather solutions below the threshold value 
(|102|) . As increases to the limit (3, power decreases, approaching both theoretical 
estimates. We also note that the lower bound ([74]) predicts the decrease of the 
power as frequency increases, approaching the threshold (|102| . 

Figure [3] considers the case (3 = 10, e = 1 in ID and 2D-lattices. We observe the 
increased accuracy of the qualitative and quantitative predictions of the variational 
lower bound ([74]) . in the ID-case. In the 2D-case an excitation threshold appears, 
i.e. there exist a minimum of the power (excitation threshold), and as frequency 
increases, the power increases, reaching a local maximum. This behaviour is in 
accordance with that described in the recent work [2D]. In the focusing saturable 
DNLS, this behaviour seems to appear in the 2D-lattice for larger values of the 
parameter /?, while in the defocusing case this behaviour seems to appear when 
only the dimension of the lattice is increased. 

The "limiting case" with respect to the size of (3, j3 = 10, clearly demonstrates 
that (|102|) is a quite sharp estimate of the threshold value on the power, for the 
existence of a breather in the focusing saturable nonlinearity. The corresponding 
breathers with large frequencies are real examples demonstrating that this estimate 
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Saturable nonlinearity. p=1 0. N=1 Saturable nonlinearity. p=1 0. N=2 
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Figure 3. Dependence of the power with respect to the frequency 
for ID and 2D lattices of the focusing DNLS with saturable nonlin- 
earity for j3 — 10 and e = 1. Power decreases approaching the lower 
bounds when £1 tends to (3. Increased accuracy of the variational 
lower bound ([74]) . for N = 1. An excitation threshold appears in 
the 2D-lattice. 



is the smallest value on the power of an existing breather. The same closeness 
appears for the variational lower bound (|74[) . 

5. DNLS equation with power nonlinearity. 

5.1. The defocusing DNLS [3 > 0: solutions ip n (t) = e- int <f> n ,n > 0. In 
this subsection we derive a threshold on the power for the existence of non-trivial 
breather solutions for the DNLS equation with power nonlinearity |T]), in the defo- 
cusing case (3 > 0. We seek breather solutions 

Mt) =e- iO *0n, fi>0, (107) 

for the DNLS equation 

iV>„ + e (A^)„-/3|^ n | 2CT ^„ = 0, 0>O. (108) 

Solutions (|107p of DNLS equation (|108|) , satisfy the equation 

- e(A d 0)„ - fl^ = -(3\K\ 2 °K, n > 0, /3 > 0. (109) 

We consider the Hamiltonian for the defocusing DNLS (| 1 08[) 

N 8 
KM = eJ2 + ^tE \M 2r7+2 , (3>0, (110) 

and the energy functional 

JV 

fibM =e^||v+0||l-n l^»l 2 ' n>0 - ( 1U ) 

We also recall that the derivative of the functional 

CM = ]T \K\ 2a+ \ 

\n\<K 
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is given by (see also [TUl Lemma 2.3]) 

{C' a [<l>]^) = (2<r + 2)Re \K\ 2 °M, for all t/j G l\7L N K )- (H2) 

|n|<iC 

Existence of a nontrivial breather solution (|107[) . will be derived by considering 
constrained minimization problems similar to those in Section 3 (see also [5T]). We 
have the following 

Theorem 5.1. A. Consider the variational problem on £ 2 (Z^) 

inf {Ha[(/)} : V[<t>] = R 2 > 0} , (113) 

for some (3 > 0. There exists a minimizer </>* €E £ 2 (Z^) for the variational problem 
H13\) and uj(R) > 0, such that f2(i?) > eXi, with both satisfying the Euler-Lagrange 
equation 

-e(A d <f>) n + P\<t> n \ 2 °<t> n = n<t> ni \n\<K, 
<j>* n = 0, \n\>K, 

with V[4>*] = R 2 - 

B. Consider the variational problem on £ 2 (Z^) 

infjfnM : Yl \M 2a+2 = M>o\, (114) 
[ \n\<K J 

for some fl > 0. Assume that 

Q > 4eA. (115) 
There exists a minimizer cf> G £ 2 (Z^) /or the variational problem ^114\ ) and j3{M) > 
0, both satisfying the Euler-Lagrange equation j!14\ ), and 53i n |<iC l ( / , n| 2<T+2 = M. 
Assume that the power of the minimizer is V[(f>] — R 2 . Then the power satisfies the 
lower bound 

"- 4 - V ' "<^:=^[# (116) 



L/3(o-+l) 

Proof: A. Relation (fTT2"J> implies that the functional 7^ : £ 2 (Z^) -> K is a 
C 1 -functional. It is bounded from below, since from (|2T))) 

AT 

> l|V+^||l > eAxi? 2 . 

The same variational arguments of Section 3 imply the existence of a minimizer 
4>* e i? 2 (Z^) of H a , and the Lagrange multiplier Q(R) > 0, such that 

JV 

(Ki^-SW^],^) = 2 £ ^(V+r,V+^) 2 + 2/3Re £ K| 2<7 #>„ 

n=i |«|<K 

-2Q ^ #$£=0, fo ra11 ^el 2 (Z^). (117) 

|n|<K 

Then setting ?/> = 0* in (|117p . and by using (|2"Tj) . we obtain that 

N 

2eAr Y ICJ 2 <26^||V+0*|| 2 + 2 / 3 ]T \^\ 2 ° +2 = 2Q £ |C| 2 , (118) 
|n|<K f=l |n|<K M<K 
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which shows that Q(R) > e\i > 0. 

B. The functional £q is bounded from below: the equivalence of norms (pi 
implies the existence of a TV-dependent constant C2, such that 

Ml < Cl\\(f>\\ 2 2a+2 , for all e ( 2 {1 N K ). (119) 

Then using (|119|) . we find that 

\n\<K 

> -fie* f x: \k\ 2 °A (121) 

> -fiCfM^. (122) 

Again the existence of the minimizer <j) and of the Lagrange multiplier \(M) 6 K 
can be obtained by the same arguments as in Section 3. Moreover by using (11121) . 
we have that 

N 

(£' n [$\-\£ R [hi) = 2 e ^(V+<£,V+V02-2mte 

v=X \n\<K 

-2(<r + l)A J2 \4>n\ 2a 4>ni^ = 0, (123) 

\n\<K 

for all tjj e e 2 (%%)- Setting ip = $ in (|T2"3"]) . we obtain 

A 

2fa[^=2e^||V+^||i-2fi £ |0„| 2 = 2(a + 1)A £ |<^+ 2 . (124) 

f=l |"|<A |n|<K 

We observe that 

£ u [4>]<AeN \M 2 -Sl J2 l^«| 2 - ( 125 ) 

|n|<K \n\<K 

Thus fn[0] < if (|1 15[) is satisfied. Note that due to the estimate (|2"2"j) . the condition 
(|115[) implies that 

fi > eAi. (126) 
Then assuming (fTT5)) . we find that A(M) < 0. We set A = ~/3, (3 > 0. Finally, we 
assume that the power of the nontrivial minimizer <fi is X^|n|<K l^™| 2 = Then, 
returning to (|124[) . and by using © which holds also in the finite dimensional lattice, 
we get 

N 

20R 2 = 2e^||V^|| 2 + 2(^ + 1)/? £ |0„| 2ct+2 

<8dV £ |^| 2 + 2(a+l)/3 ]T |0„| 2ct + 2 
|n|<K |n|<K 

<r+l 

<8dV Y, l^| 2 + 2(a+l)/3l 2 p 

\n\<K \\n\<K 

<8eNR 2 + 2{a + l)(3R 2a+2 . (127) 
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Under condition (|115j) . inequality (|127| implies the lower bound (I116j) . o 

A threshold for the power could be also derived from the case A. of Theorem 15. II 
working exactly as for the derivation of that in (|127[) . we find from (|118p , that 

Sli? 2 < AeNR 2 + (3R 2a+2 . (128) 

Thus, in the case of Theorem 15.11 A. and under the hypothesis that the Lagrange 
multiplier Sl(i?) of the case A is taking values SI > AeN, we find from (|128|) that 

" <R 2 :=V 2 W]. (129) 

Comparing with (|116| , it readily follows that 

Vi[4>]<V 2 [4>*]. (130) 

5.2. Numerical study for the defocusing DNLS with power nonlinearity. 

Similarly to the results of Section 3, for the saturable DNLS, it seems interesting 
to test the behaviour of the lower bounds (|116|) and (|129p . as thresholds for the 
existence of breather solutions for the defocusing DNLS (|108p . Theorem 15.11 A 
implies that, for a given (3 > 0, the existence of a frequency SI > eAi (as a Lagrange 
multiplier), and of a nontrivial minimizer (jf of the Hamiltonian (|110|) . such that the 
corresponding breather solution ip n (t) = e~ l!lt 0n has a power satisfying the lower 
bound (|129p . in the case where f2 is assumed to be such that SI > AeN . On the other 
hand, Theorem l5.1I B implies, for a given SI > 4eN, the existence of /3 > and of a 
nontrivial minimizer <p of the energy functional (jllip . such that the corresponding 
breather solution ip n (t) — e~ int <fi n , has a power satisfying the lower bound (|116p . 

The numerical study is for ID and 2D-lattices. For the case N = 1 we consider 
e = 0.25, and for the case N = 2 we consider e = 0.15. We study first values 
of a satisfying a > j^. In Figure [?J the study refers to the cases a = 1, N = 2, 
and a = 2, N = 1 respectively. We observe first that the numerical power of the 
solutions fulfils both lower bounds (|129|> . (|116p . Moreover we observe that they can 
be considered also as thresholds on the power of periodic solutions with frequencies 
SI > AeN. This fact is revealed by the case a = 10, N = 1 for which the lower bound 
(|116p is proved to be a quite sharp estimate of the power for large frequencies. The 
same satisfactory accuracy of the theoretical estimates (|116p and (|129p is observed 
also in case ofcr = 10,iV = l considered in Figure [5] We note that the phonon band 
of the defocusing DNLS equation extends to the interval [0,4eiV]. Then breathers 
frequencies must lie in the intervals SI > 4eN, or SI < 0. It is the former case 
which we consider in this paragraph. The numerical studies for the case a < 2/N 
where the excitation threshold [21] do no exist, seem to fully justify the argument 
that the lower bounds (|129p . (|116p can be used as thresholds on the power for 
the existence of breather solutions. The results for a < 2/N are demonstrated in 
Figures [6] and [7l Figure [6] considers the case a = 0.1 for N = 1 and N = 2 and 
(|129p , (|116p are clearly sharp estimates of the smallest value of the power a breather 
solution can have. This is justified by the numerical power of breather solutions 
with small frequencies. The numerical power for the case a = 1, N = 1 also fulfils 
the theoretical estimates. 

We remark that the approach of minimizing the linear energy appeared to be 
useful also in the defocusing DNLS with power nonlinearity since the lower bound 
(|116p provides in general better quantitative predictions of the numerical power if 
compared with (|129p . 
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Defocusing power nonlinearity.a=1 . N=2 



Defocusing power nonlinearity. c=2. N-1 




Figure 4. Numerical power for the defocusing DNLS with power 
nonlinearity and (a) a = 1, N = 2, e = 0.15 (b) a = 2, N = 1, e = 
0.25. The lower bounds (|116p and (|129[) . dot-dashed and dashed 
lines respectively, serve as thresholds for the existence of nontrivial 
breather solutions. 



Defocusing power nonlinearity. g=10. N=1 




Figure 5. Numerical power for the defocusing DNLS with power 
nonlinearity and a = 10, N = 1, e = 0.25. Both lower bounds ()116j) 
and (|129[) . dot-dashed and dashed lines respectively, are fulfilled as 
threshold values. 



5.3. Focusing DNLS with power nonlinearity ((3 < 0): solutions 4> n (t) — 
e~ iat (j) n , fi < 0. We conclude with a numerical study of the focusing DNLS with 
power nonlinearity. Thus we shall consider solutions with frequencies 51 < (the 
latter case of the phonon band condition). As in Section 3, we use for convenience, 
the notation (3 = —A, A > 0, and £1 = — ui, uj > 0. 

Instead of the threshold estimate [TDl (2.38), pg. 126] derived by the fixed point 
argument, for breather solutions 



ip n (t) = e^Vn, ^>0, 



of the DNLS 



in infinite lattices, 



mm 



l 



A(2cr + 1) 



(131) 



(132) 
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Defocusing power nonlinearity. a=0.1 . N=1 Defocusing power nonlinearity. a=0.1. N=2 




Figure 6. Numerical power for the defocusing DNLS with power 
nonlinearity case for (a) cr = 0.1, N = 1, e = 0.25 (b) a — 0.1, 
N = 2, e = 0.15. The lower bounds (fTTB|> and (fT2"9"j) . dot-dashed 
and dashed lines respectively, arc fulfilled as threshold values. 




Figure 7. Numerical power for the defocusing DNLS with power 
nonlinearity case a = 1, N = 1, e = 0.25. The lower bounds (|116p 
and (|129p , dot-dashed and dashed lines respectively, are fulfilled as 
threshold values. 



we may derive a possibly improved one: by using (|92p (which holds also for an 
infinite lattice) instead of [TOJ, (2.44), pg. 127], we find that 



Pu ■ = 



A(2a+ 1)_ 



= 4E Z 



(133) 



Although (|133p . holds also for the case of the Dirichlet boundary conditions, we 
may derive a threshold value taking into account the finite dimensionality of the 
problem: for solutions (|13ip . we consider the operator 7o : ^ 2 (Z|£) — > i? 2 (Z^), 
defined as 

{Tu<j))\n\<K = -z{&d<l>)\n\<K + U<t>n- (134) 
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Focusing power nonlinearity. a=1 . N-1 Focusing power nonlinearity.a-10. N=1 




Figure 8. Numerical power for the focusing DNLS with power 
nonlinearity and (a) a = 1, N = 1, e = 1 (a < -w-no excitation 
threshold ), (b) a = 10, N = 1, e = 1 (a > -^-excitation threshold). 
The inset in (b) shows a magnification of the region where the 
power reaches its minimum value. Dashed line corresponds to lower 
bound (fT37l) 



Clearly 



Si := eXi+uj > 0, 



and by using (|20|) . we get that 



N 



■u\\<i>\\l > tfi||0||a for all <jy E f(Z%) 



Then, by using (f93|) . we may derive the threshold value 

eAi + uj 



(135) 



(136) 



(137) 



_A(2a+ 1) 

Throughout this study we have fixed e = 1. A first result is that the power of 
the solutions fulfils (|133[) , that is, it is higher than 



_A(2er + 1)_ 

We mention first, that regarding the numerical study of the threshold (|137p . as in the 
case of the saturable nonlinearity, we have not observed remarkable improvement or 
differences in comparison with (| 1 33[) . Although the threshold (| 1 33[) is independent 
of the dimension, it is interesting to compare this with the results of [21], related to 
the conditions on existence of excitation threshold which depends on the dimension 
and the nonlinearity exponent a. In Figure [51 we present the numerical power 
(solid curve) for the cases a = 1 and a = 10, and N = 1. The values a — 1 and 
N = 1 satisfy a < 2/N, and the power fulfils the threshold (|133[) for existence. 
The case a = 10 and N = 1 satisfies a > 2/N, and as predicted in [5T] and [B], 
the power approaches a minimum (the excitation threshold) . The numerical power 
still fulfils the threshold (|133p . In the latter case, the corresponding breathers with 
large frequencies provide real examples demonstrating that Q133[) is a quite sharp 
estimate of the smallest value of the power a breather can have when a > 2/N. 
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Focusing power nonlinearity. a=0.1 . N=1 Focusing power nonlinearity.a=1 . N=2 




Figure 9. Numerical power for the focusing DNLS with power 
nonlinearity and (a) a = 0.1, N = 1, e = 1 (cr < ^-no excitation 
threshold ), (b) a = 1, N = 2, e = 1 (cr = -^-excitation threshold). 
The inset in (b) shows a magnification of the region where the 
power reaches its minimum value. Dashed line corresponds to lower 
bound (fm]) 



Focusing power nonlinearity. a=0. 1 . N=2 Focusing power nonlinearity. a=10. N=2 




Figure 10. Numerical power for the focusing DNLS with power 
nonlinearity and (a) a = 0.1, N = 2, e = 1 (cr < ^--no excitation 
threshold ), (b) a — 10, N = 2, e = 1 (cr > -^-excitation threshold). 
The inset in (b) shows a magnification of the region where the 
power reaches its minimum value. Dashed line corresponds to lower 
bound (fj~37l) 



Note that the dependence of this excitation threshold with respect to cr has been 
numerically calculated in [13]. 

In figure [5] we present the results for the cases a = 0.1, N = 1 and cr = 1, N = 2. 
The first case satisfies a < 2/N and the second satisfies cr > 2/N, in its critical 
value cr = jj. The threshold (|133[) is fulfilled in both cases. We observe that in the 
case a = 0.1, N = 1, the threshold (|133[) is a sharp estimate of the power observed 
for small values of frequency. These breathers serve this time as real examples 
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demonstrating that (I133j) predicts the smallest value of a breather also in the case 
(T < 2/N. In the second case we observe the appearance of the excitation threshold. 

To show the global character of the estimate (| 133[) as a threshold on the existence 
of breather solutions, we consider in Figure [TOl the "limiting" cases with respect to 
the size of the nonlinearity exponent, a = 0.1 and a = 10, this time for N = 2. 
The first case is again an example for a < 2/N and the second for a > 2/N. The 
threshold (|133p is again fulfilled in both cases. We observe that in both cases the 
threshold (|133[) is still a quite sharp estimate of the smallest value of the power 
of an existing breather as the comparison with the power of breathers with small 
(large) values of frequency demonstrates. 

It is worth noticing the existence of a maximum in the power for a = 1, N = 2, 
as was predicted in |12j . 
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